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ABSTRACT 

In this report, a robust error control coding scheme is presented. The 

scheme is a cascaded FEC scheme supported by parity retransmissions for further 

error correction in the erroneous data words. The error performance and 

throughput efficiency of the scheme are analyzed. Two specific schemes are 

proposed for NASA near earth satellite communications. We show that both 

schemes provide high reliability and throughput efficiency even for high channel 

*2 

bit -error rates in the range of 10 . The schemes are suitable for high data 


rate file transfer. 


PERFORMANCE ANALYSIS OF A HYBRID ARQ ERROR CONTROL SCHEME 
FOR NEAR EARTH SATELLITE COMMUNICATIONS 


1. Introduction 

In an earlier technical report [1], we proposed a hybrid ARQ error control 
scheme for NASA's near earth satellite communications. The scheme was designed 
to provide very low error probability and high throughput for high data rate 
file transfer. In this report, we analyze this error control scheme and show 
that it indeed performs very well in terms of error probability and throughput 
efficiency, even for very high bit-error-rates (BER) , say in the range 10 

In order to make this report self-contained, we re -describe the proposed 
scheme here. The proposed error control scheme is a FEC (forward-error- 
correction) scheme supported by parity retransmissions for further error correc- 
tion [2 , 3 ,4] , . known as a hybrid ARQ scheme. The FEC scheme is achieved by 
cascading two linear codes, the outer and inner codes. The inner code is either 
a binary block code or a binary convolutional code. The outer code is a block 
code with symbols from the Galois field GF(2 m ), say a Reed-Solomon (RS) code. 

The retransmission is designed to provide additional parity-check symbols for 
further error correction. These additional parity-check symbols are used 
without decreasing the rate of the overall cascaded code. 

The report is organized in the following manner: the FEC scheme is des- 

cribed in Section 2, the retransmission strategy is presented in Section 3, 
performance analysis of the proposed scheme is given in Sections 4 and 5, 
specific schemes using NASA standard RS outer code are given in Section 6, and 
conclusion remarks are given in Section 7. Computation results on error proba- 
bility and throughput performance are included in tables and figures. 

2 . The FEC System 

The FEC scheme is a cascaded coding scheme [5,6]. In this study, the inner 
code, denoted C^, is chosen as a binary (n^,k^) linear block code with minimum 




» 


distance d^. is designed to correct t^ or fewer errors and simultaneously 

detect A^ (with A^ > t^) or fewer errors where t^+A^+1 < d^. The outer code, 

denoted C Q , is obtained by interleaving a maximum-distance-separable (^,1^) 

2 

linear code C 2 over GF(2 ) with minimum distance d^ and rate greater than 1/2. 

Let m^ be the interleaving degree (or depth) . Then the outer code of the scheme 
l 

is an (m^n 2 .m.^) code with symbols from GF(2 ). Note that - ^-l^+l [2]. 

We call C£ the base outer code which is designed for correcting symbol errors 
and erasures. In our proposed scheme, we assume that the code parameters 
satisfy the following conditions: 

k^ - m^X 

( 1 ) 

n 2 - m 2 (n 2 -k 2 ) - m 2 (d 2 *l) , 

with m^ > 1 and m 2 > 2. 

There is a third code which is used for parity retransmissions. C r is a 

half-rate (2d 2 *2, d 2 *l) code obtained by shortening the base outer code C^. 

is also a maximum- distance -separable code with the same minimum distance d 2 as 

C,,. Since C is a shortened code obtained from C„, C„ and C can be encoded and 
2 r 2 2 r 

decoded by the same circuits. A very important property of C r is its invertible 

structure which will be used in data recovery during the parity retransmissions. 

— 2 — 

Let u be a sequence of d 2 *l information symbols from GF(2 ). Let R r (u) denote 

the sequence of d 2 *l parity-check symbols formed based on the information 

sequence u and the half-rate code C^,. Then the 2(d 2 *l) -symbol word (u,R r (u)) is 

a codeword in C^. There is a one-to-one correspondence between u and R^(u). 

Consequently, knowing only the d 2 -l parity-check symbols of a codeword in C^, 

the corresponding d 2 *l information symbols can be uniquely detrmined by an 

inversion operation on the d 2 -l P ar ity-check symbols [see Appendix A]. If C is 

a shortened cyclic code, the inversion from R^(u) to u can be achieved by a 

feedback shift register [2]. Let v^ be a codeword in C 2> Since n 2 - m 2 (d 2 -l), 
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we can divide v into m 0 subsequences, v. 1( v. 0 v. ; each consists of d„-l 

1 ^ ii 1/ i,®2 l 

symbols. For 1 iS j < n^, let R r (v^) be the sequence of c^-l parity-check sym- 
bols formed based on and C r> Clearly (v^ ,R r (v^ ) ) a codeword in C^. 

Let 


R( V - <\<’ u ).* r (v l2 ) V v i.m 2 » 


( 2 ) 


Then it can be shown that R(v^) is also a codeword in C 2 [see Appendix B] . In 
fact, is a codeword in C 2 if and only if K(y^) is a codeword in C^. This 
property will be used in our proposed error control scheme. For convenience, we 
call R(v^) the parity word of v^. 


Encoding 

A message consists of a string of ^ x k 2 information bits. This string is 
divided into k 2 segments , each segment consists of k^ information bits. Each 
segment is further divided into m^ i- bit bytes. Each i-bit byte is regarded as 

a 

a symbol in GF(2 ). The encoding operation consists of two stages as shown in 
Figure 1. For each input message of k^k 2 kits, the output is an n^n 2 -bit code- 
word in the cascaded code C. A codeword in two-dimensional format is shown in 
Figure 2. The transmission is done column by column and from left to right. At 
the first stage of encoding, each k^-bit segment is encoded into an n^-bit code- 
word in the inner code C^, which Is called a frame. At the same time, the m^ 
i-bit bytes are multiplexed into C^-code encoders to form parity-check sym- 
bols for m^ codewords in C 2 « As soon as the k 2 segments of a message have been 
shifted into the overall encoder, k 2 frames have been formed and transmitted. 
Also all the parity- check symbols of codewords in C 2 have been formed and are 
in the registers of the m^ Opcode encoders. Then these parity-check symbols 
(m^(n 2 -k 2 ) of them) are multiplexed and shifted into the inner code encoder to 
form n 2 *k 2 more frames (they are parity frames). These n 2 -k 2 parity frames and 
the k 2 data frames formed at the first stage together form a complete codeword 
array in the cascaded code C. 


- 3 - 


> 


There is another part of the encoder. This part is for retransmission (if 
needed). It consists of an encoder for the half-rate code C r> and a buffer. 

Let v i' v 2 ,, * ,,v m be codewords in C ^ which are formed by the upper part of the 

overall encoder. The function of C^-encoder is to form the m^ parity codeword 
in C 2 , 

R(v ) , R(v ) , .... R(v ) 
l z 

corresponding to the codewords, v ]^ » v 2 » ’ * • * v m > where R(v^) is given by (2). 

These m^ parity words are temporarily stored in a buffer for possible 
re tr ansmi s s ion . 

There is another encoding arrangement. We can use a single C ^- encoder to 
form m^ codewords in (the first rows in Figure 2) and store them (in an 
array form) in a buffer. Then encode the columns into n^ frames, and 
transmit them column by column. This encoding arrangement requires more buffer 
store . 

Corresponding to each message of k^ x k^ information bits, the output of 
the overall encoder shown in Figure 1 is a string of n^ frames. These n ^ frames 
are said to form a data-block. Consider the data-block shown in Figure 2. The 
top k^ - m^i rows of the array is actually regarded as m^ i-bit byte rows. Each 
of these i-bit byte rows is a codeword in aTU * called a data-section of the 
code array. The m^ data-sections form a subarray which is called a data-segment 
array. Each column of a data- segment array is a data- segment . There are k£ 
message segments and n^-k ^ parity segments. Each data section is further 

divided into m^ subsections, each subsection consists of s y m ^ > °^ s ( or i-bit 

i - - - 

bytes) from GF(2 ). The m parity words, R(v ) , R(v 9 ) , ..., R(v ) form a 

X 4 ^1 

parity- segment array in the buffer. Each column of this array will be called a 
parity-segment. Each row is called a parity-section, and consists of m 2 
subsections . 
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Decoding of a Data Block 


The decoding of a data-block is basically the same as the one described in 
our earlier technical report on, "A Cascaded Coding Scheme for Error Control" 

[5]* It consists of two stages. The first stage of decoding. Depending on the 
number of errors in a received frame, the inner code decoder performs one of the 
three following operations: error -correct ion , erasure and leave- it-alone (LIA) 

operations. When a frame in a data block is received, its syndrome is computed 
based on the inner code C^. If the syndrome corresponds to an error pattern e 
of t^ or fewer errors, error correction is performed by adding e to the received 
frame. The n^-k^ parity bits are removed from the decoded frame, and the 
decoded m^-byte data- segment is stored in a receiver buffer for the second stage 
of decoding. A successfully decoded data- segment is called a decoded segment 
with no mark . Note that the decoded segment is error- free , if the number of 
transmission errors in a received frame is t^ or less. If the number of trans- 
mission errors in a received frame is more than the errors may result in a 
syndrome which corresponds to a correctable error pattern with t^ or fewer 
errors. In this case, the decoding will be successful, but the decoded frame 
(or segment) contains undetected errors. If an uncorrectable error pattern is 
detected in a received frame, the inner code deocder will perform one of the 
following two operations based on a certain criterion: 

1. Erasure Operation - The erroneous segment is regarded being erased. In 
fact this segment is not really removed from the buffer, it is still 
stored there for later use. This segment is called an erased segment. 
Each i-bit byte of an erased segment is regarded as an erasure for the 
outer code decoding. 

2. Leave- it- alone (LIA) Operation - The erroneous segment is stored in the 
receiver buffer with a mark , We call such a segment a marked segment. 
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Thus, after frames of a received block have been processed, the receiver 
buffer may contain three types of segments: decoded segments without marks, 

erroneous segments with marks , and erased segments . 

The above inner code decoding consists of three operations: error- 

correction, erasure and LIA operations. An inner code decoding which performs 
only the error -correction and erasure operations is called an erasure- only 
decoding . On the other hand, an inner code decoding which performs only the 
error-correction and LIA operations is called a LIA-only decoding. 

When frames in a received data-block have been processed, the decoder 

buffer contains a decoded data-segment array with m^ i-bit byte rows. Each of 

these i-bit byte rows is regarded as a received codeword from which may 

contain erroneous symbol s (marked or unmarked) and erasures. The code C ^ and 

its decoder are designed to correct the combinations of symbol erasures and 

symbol errors. Maximum-distance-separable codes (or Reed-Solomon codes) with 
1 

symbol from GF(2 ) are most effective in correcting symbol erasures and errors. 

At the second stage of decoding, the C^- decoder attempts to decode the rows 
of the data-segment array. Let i and h be the numbers of erased and marked 
segments respectively. The receiver stops the decoding process and requests a 
retransmission for the erroneous data-block if either of the following two 
events occurs: 

1. the number i is greater than a certain pre-designed erasure threshold 

T with T < d 0 -l. 
es es 2 

2. the number h is greater than a certain pre-designed threshold T^(i) 
with T ei (i) < L(d 2 -l-i)/2j for a given i. 

If none of the above two events occurs, the -decoder starts the 
error -correction operation on the erroneous sections (or rows) of the data 
segment array, one at a time (they can be processed at the same time if we use 
m^ C 2 -decoders) . The i symbol erasures and the symbol errors with or without 
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marks in each section are corrected based on the code C 2< Let t 2 (i) be the 
error -correct ion threshold for a given i where 

t 2 (i) < L(d 2 -l-i)/2j . (3) 

If the syndrome of a section in the data- segment array corresponds to an error 
pattern of i erasures and t 2 (i) or fewer symbol errors, error -correction is 
performed. The values of the erased symbols, and the values and the locations 
of symbol errors are determined based on a certain algorithm. If more than 
t 2 (i) symbol errors are detected, then the receiver stops the decoding process 
and requests a retransmission for the erroneous data-block. If all the m^ sec- 
tions of a data segment array are successfully decoded, then the k 2 decoded 
data- segments are either delivered to the user or saved in the buffer until they 
are ready to be passed to the user. 


3 . Retransmission Scheme 

When the receiver fails to decode a data-block v, it saves the erroneous 

data-segment array of v in a buffer and requests a retransmission for v. The 

retransmission is not v itself but a parity-block P(v) corresponding to v. The 

parity-block P(v) is formed based on v. Let v. ,v_ v be the m, sections of 

I L m^ 1 

the data segment array of v. Recall that each section v^ is a codeword in C 2 » 
For each section v^, the encoder has already constructed a corresponding parity 
codeword 


R(v.) - (R (v._), R (v._), ..., R (v. )) 

i r ir r v t2' ’ ’ r v i,m 2 

in C 2 where R^(v„) is the parity-check part formed based on the j-th subsection 

v. of v. and the half- rate (2d. -2, d„-l) code C (i.e., (v..,R (v. .)) is a 
1 j i 2 2 r ' 1 j r v i j ' 

codeword in C^) . The m^ parity codewords, R(v^) ,R(v 2 ) , . . . ,R(v m ) are stored as 

a m^ x n 2 segment array in the transmitter buffer, which is called a parity- 
segment array. When the transmitter receives a request for a retransmission for 
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data -block v, the inner code encoder encodes each segment of the parity- segment 
array into a frame (a codeword in C^) . Hence a parity-block P(v) is formed, 
which is also a codeword in the cascaded code C. The parity-block P(v) is then 
transmitted to the receiver. 

When a parity-block P(v) is received, the receiver starts to decode it. 

The decoding of P(v) is the same as the decoding of a data-block v. If the 
decoding of P(v) is successful, inversion is then performed on the first k£ 
decoded segments of the parity- segment array. This inversion gives the k£ 
data- segments of v. These decoded data-segments are then delivered to the user 
or saved in the receiver buffer until they are ready for delivery. At this 
time, the erroneous data-segment array which is stored in the receiver buffer is 
discarded. 


If the C2*decoder fails to decode the parity-block P(v) , then the parity- 
segment array of P(v) and the data-segment array of v (which is stored in the 
receiver buffer) together are used for error correction based on the half- rate 


code C^. The receiver puts v„ and R r (v„) together to form a subsection pair, 

(v. . ,R (v..)). Then the C -decoder decodes (v. . ,R (v..)) into an estimate v * . 
ij r ij r ij r ' ij 7 ' ij 

for v £j • After x m2 such decodings, the receiver contains the following 
estimated data-segments array: 



Then the receiver checks whether each v^, for 1 < i < m^, is a codeword in C2. 

^ 'Jc ^ 

Note that this time C„ is used only for error detection. If all v, ,v„,...,v 

Z 1 Z m. 
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are codewords In C^, then the decoding is successful and the k 2 estimated 
data- segments are accepted by the receiver. If any v^ is not a codeword in , 
the receiver discards the erroneous data- segment array of v (stored in the 
buffer) and save the parity- segment arragy of P(v) for later use. At the same 
time the receiver requests a second retransmission for v. The second retrans- 
mission is the data-block v itself. When v is received, it is decoded as 
before. If the (^-decoder fails to decode v, then the data-segment array of v 
and the parity- segment array of P(v) (stored in the buffer) together are used 
for error correction based on C^. If the correction process is not successful, 
then the receiver discards the parity- segment array of P(v) and saves the data- 
segment array of v. At the same time, the receiver reqeusts a third retransmis- 
sion for v. The third retransmission for v is the parity block P(v) . When P(v) 
is received, the receiver starts the decoding process again. Therefore, the 
retransmissions are alternate repetitions of the parity block P(v) and the 
data-block v as shown in Figure 3. The retransmissions continue until the 
message in v is finally recovered by the receiver. 

The major advantage of this error control scheme is that extra parity 
symbols for error correction are transmitted only when they are needed. These 
extra parity symbols are used without decreasing the rate of the cascaded code 
C. If the half-rate code C is powerful enough, at most one retransmission is 
needed to recover a message. When the channel is not very noisy, the error 
correcting capability of the cascaded code C should be able to recover the 
message in its first transmission. In this situation, the system throughput is 
equal to the rate of C which is R^R 2 . wtlen the channel is noisy, the first 
retransmission provides us the parity symbols of for extra error correction 
capability. Since is used for correcting errors only in a subsection of a 
codeword in C^, and since C has the same error correcting capability as C 2 > 
errors in the entire data-segment array should be corrected by C^. In this 
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situation, the throughput of the system should be If the noisy situa- 

tion is rare , then the proposed error control schemes provides maximum 
throughput most of the time. 

Since is a shortened code obtained from C^, the decoder for C ^ can be 
used for decoding C . Therefore, only decoders for and C ^ are needed. Since 
the inner code is binary and shorter, its decoder is much simpler than the 
decoder for C^. 

In our earlier report [5], we showed that a cascaded coding scheme provides 
extremely high reliability. We expect the proposed scheme in this report will 
also provide extremely high reliability. Analysis of the scheme will be given 
in sections 4 and 5 . 

A special case for the above error control scheme is that n^ - *• i. In 

this case, no inner code is used, the outer code is simply C ^ which is used for 

both error correction and detection. The code C is used for error correction 

r 

only. 

Retransmissions can be carried out in any of the three modes: the stop- 

and-wait, the go-back-N and the selective -repeat [2]. Selective -repeat scheme 
is the most efficient retransmission scheme and provides the highest throughput 
efficiency. For high data rate file transfer over satellite links, only 
selective -repeat ARQ scheme provides satisfactory throughput for high channel 
bit-error rate. Of course, selective -repeat ARQ scheme is more complicated to 
implement than the other two ARQ schemes. Various selective -repeat ARQ proto- 
cols have been proposed and studied [2,7-10]. A practical and efficient 
selective-repeat ARQ protocol, particularly designed for satellite communica- 
tions, is the Yu-Lin's selective -repeat ARQ protocol [2,7]. The basic Yu-Lin's 
protocol requires a storage buffer at the receiver which is capable of storing N 
blocks (data or parity) that can be transmitted during a round-trip delay 
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period. This basic Yu-Lin's protocol has been proved to provide very good 
throughput performance for satellite links. Better throughput performance can 
be achieved by doubling the size of buffer store [7]. 


4. Performance Analysis - A Special Case 

In this section, we analyze a special case of the proposed error control 
scheme for which n^ - k^ - i and m^ - 1. That is, in this special, there is no 
inner code and the base outer code C£ is not interleaved. Hence every data 
(or parity) block consists of a single codeword from C^. Since there is no 
inner code, there is no erasure operation. The outer code C£ is used to correct 
up to t£ ^ symbol errors. The half-rate code is used to correct 

up to t^ - L(d 2 *l)/ 2 j symbol errors in each data-parity subsection pair when 
retransmission occurs. To evaluate the overall error performance of the scheme, 
we need to know the error performance of the C 2 *decoder and C^-decoder. 


Error Performance of decoder 

Let € be the channel bit error rate. Let P and P denote the probabil- 

er es r 

ities of an incorrect decoding and a decoding failure respectively for a 


received data (or parity) block by the C 2 'decoder. Then 

vJ 


and 


where 


P + P 
er es 


- 1 LiM'KP 

j-t.+l u J 


n 2 min{t2,n 2 -w) 


er 


-5 


w— d„ 


n„ 


h-0 


n 2 -w 


U J 


j-w+h-t 2 


P e - 1 - (1-0 . 


w+h-d 


P(w,h,j) - (p ] 


2 „ - 
U-P e ] 


5 

i-0 


(4) 


(5) 

( 6 ) 

(7) 


The derivations of (4) and (5) are given in [5,6]. 
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Let P denote the right-hand side of (5). Then P^„ is an upperbound on 
er er 

the probability of a decoding error of C 2 -decoder. In the first specific scheme 

proposed for NASA near earth satellite communications, the (256,224) RS code 
o 

over GF(2 ) is chosen as the C 2 code. The code is capable of correcting a 
maximum of 16 symbol errors. For several values of t. < 16, P and P +P are 

Z w L w JL CD 

given in Table 1 and Figures 4 and 5 . 


Error Performance of C -decoder 
r 

Let v and R(v) be a received data-block and its corresponding parity-block 
respectively. For 1 < j < m 2 , let v^ and R^(Vj) be the subsections of v and 
R(v) respectively. If v. and R r ( v j) are error-free, the word (v^ , R r ( v j)) is a 
codeword in C^. In a retransmission, when the C ^- decoder fails to decode R(v) 
(or v) , the C^-decoder starts to decode the m 2 subsection pairs, (Vj,R^(v^)) for 

>n> 

1 < 1 < m„. The C -decoder decodes (v. ,R (v. ) ) into an estimate v. if the syn- 
J 2 r j r j j J 

drome of (v^.R^Vj)) corresponds to [(d 2 -l)/2j or fewer errors. If the 

C^-decoding is successful for every subsection pair and the decoded block 

-* -* -* -* 

v - (v lt v 2 v_ ) 


mo 


is a codeword in C 2> then v is delivered to the user. In this case we say that 

the C r *decoding is successful for the data-block v. Otherwise, a decoding 

failure is declared, and another retransmission is requested. To achieve high 

reliability, the base outer code C 2 should be a rather powerful code (like the 

8 

NASA standard (255,223) RS code over GF(2 )). As a result, is extremely 

powerful because it is used only to correct errors in a subsection of v. We can 

imagine that the probability of a decoding failure of C^-decoder is very small 

even for a high bit-error rate. Therefore, for practical applications, we may 

assume that at most one retransmission is needed. 

Assume that a data-block v and its parity-block P(v) are sent and the 

( r) ( r) 

received pair is decoded by the C -decoder. Let P and P' be the 
r r er dr 
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probabilities of a decoding error and a decoding failure of the C^-decoder 


respectively for a subsection pair (v^.R^/v^)) from v and P(v) 


n 


P (r) + P*J* 

er df 


j-t +1 
J r 


r "\ 

n 

r 


U J 


- i J r. - 1 v J 


i?.rH 


n 


p (r) < ? (r) 
er er 




r f n ' 


w-d„ 


w 


min{t r> n r -w} 


h-0 


n -w 
r 


I 0 

i-w+h-t 


Then 


P(w,h, j) 


( 8 ) 


( 9 ) 


where n 2 - 2(d 2 *l) and t^ - [(d 2 -l)/2j . Let P^^^^ and P^f 1 ^ be the probabil- 
itles of a decoding error and a decoding failure of the C^- decoder respectively 


for the data block v. Then 


p(an) + (all) p(r> . 

er df L er df J 


df - * ■ i* ~ i • (i°) 

Note that, for any combination of symbol errors introduced by the 
C^-decoder in the first n^-l su ^sections of v, there is exactly one symbol error 
pattern introduced by the C^-decoder in the last subsection of v (which consists 
of all the check symbols) such that the resulting block v is a codeword in C^. 
The probability that a given nonzero error pattern is introduced by the 
(^-decoder in a specific subsection of a data-block is upper bounded by 


d 2 -2 


p(max) 
ep 


5 

h-0 


j-L(d,+l)/2j 


t' ! ) [-.]>■*.] 


V2 ' h 


(ID 


[see Appendix C for proof] . Then we have 

V 1 




,-i-l 


(12) 


i-1 


Overall Error Performance and Throughput Efficiency 

As we pointed out earlier, since is very powerful, at most one retrans- 
mission is needed. Therefore, the performance with at most one retransmission 
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will give us a good indication of the overall performance of the system. Let 

( 2 ) ( 2 ) 

and be the probabilities of a decoding error and a decoding failure 

for the system. Let denote the right-hand side of (12). Then we obtain 

the following upper bounds on the error performance of the system: 

P (2) < P + P P + P (all) < p + (P +P )P + P (all) a P (2) , (13) 

er er es er er er er es er er - er 


P<’> < P<» 11 > S P<» ll > ♦ P<f » A P< 2 > 
df df er df - df 


.( 2 ) 


Define 17 as follows: 


(14) 


.(2) A 


V 1 +P es> 


.( 2 ) 


(15) 


Then t) is a measure of the system throughput efficiency. Obviously 77 
lowerbounded by 

( 2 ) n 2 - d 2 +1 


( 2 ) 


is 


n„(l+P +P ) ' 
2 V er es 7 


(16) 


5. Performance Analysis - General Case with Erasure-onlv Decoding 

In this section, we analyze the performance of the proposed scheme in which 
the inner code decoder performs only the error-correction and erasure operations 
(no LIA operation) . 

Error Performance of the Combined and C ^ Decoders 

In each transmission (or retransmission) , a data (or parity) block is first 

decoded by the C^-decoder and then by the C 2 *decoder. The C^-decoder may result 

in erased segments, each erased segment contains m^ erasure symbols. The 

C 2 - decoder then attempts to correct the erasures and symbol errors in each 

section of the received data (or parity) segment array. Again let P gr and P gs 

denote the probabilities of a decoding error and a decoding failure for a block 

respectively by the combined and C 2 decodings. An upper bound on P and an 

expression for the sum P +P have been derived and can be found in one of our 

er es 

earlier technical reports [5]. 
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Error Performance of the C^- decoder 


Let v and P(v) be a received data-block and its corresponding parity block 

respectively. When the C 2 - decoder fail to decode v and P(v) in two separate 

transmissions, the C -decoder starts to decode. At this time, the receiver 

r 

buffer contains two segment arrays , one obtained from v and the other obtained 

from P(v) . Both segments are obtained after the inner code decoding, and hence 

both contain erased segments. For 1 < i s m^ and 1 < j < let v^. and 

R r (v^j) be the j-th subsections of the i-th sections of the data-segment array 

and its corresponding parity- segment array respectively. Let s^ and sj be the 

number of erasures in v^. and R^Cv^) respectively. For each subsection pair 
— — 

(v^ . , R r (v^.)), the (^-decoder decodes it into an estimate v^., if its syndrome 


correspond to Sj+sj symbol erasures and [^-l-s^ -sj )/2j or fewer symbol errors. 


-* -* _* 


If the C^-decoder successfully decodes every subsection pair, and v^ - ^ v il’ v i2 


-* 


. . . ,v. ) is a codeword in C„ for every i with 1 ^ i < m n . Then the decoded 

1.^2 2 1 


words, v i> v 2 >‘*’> v m , with their parity symbols removed are delivered to the 

user, and the (^-decoding is said to be successful. Otherwise, a decoding 
failure is declared and another retransmission is requested. 

Now we analyze the error performance of the C - decoder. Assume that a 
data-block and its corresponding parity-block are sent and the received pair is 
decoded by the (^-decoder. For 1 £ m^, let P^^ (u) and P^Jp(u) be the prob- 

abilities of a decoding error and a decoding failure respectively for a subsec- 
tion pair of the u-th section of a data-parity block pair by the C^,- decoder. 

( 1 ) l 

Let denote the probability of a frame erasure. For any element a in GF(2 ) 

and 1 < u < m^, let p^Cu.a) be the probability that a segment is not erased and 

the u-th error symbol of the decoded segment is a. For 1 < u < m^, let p g (u) be 

the probability that the u-th symbol of a decoded segment is erroneous. Then 
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( 17 ) 


P e (u) “ 


P (u,a) 


aeGF (2 ) 

Qf^O 


A procedure for evaluating p e (u,a) is given in [6], Using the same argument as 
in [6] , we can show that 


r )( “) + P d£ >(u) - } Li Kl’] 1 t ' r - 1 

i-0 L1 J j - [n - l- i)/2j+l'- J 


■ [p e (u)] J [l-P<J>.5 e (u,] n ' 1 3 


) " r [p^l 1 [i-p' 1 )! 

L • L es J L es . 

• j 


n -l 
r 


n -i mini [(n -l-i)/2j ,n -i-w) 

1 -- 1 


i?" «■£?«• } [;i }■ f":’ 1 ] e 

i-0 ^ ' w-d„-i'- w ' h-0 


where n r - 2(d 2 *l) and 


r 2 f W ] P(u * i w h *J ) 

h J j-w+h-L(n r -l-i)/2j uJ 


- n\]ir -ii+w+h-d- r -,n 

P(u,i.w ( h,j) - ^ s ; j [p e (u)j 2 |^p e (u, 0 )j r 


in -i-w-h 


[-"T" 5 k.-'.] 


j+d 2 -i-w 


and 7 is a primitive element in GF( 2 i ). Let p^ a11 ^ and be the 

er df 

probabilities of a decoding error and a decoding failure for the data block v. 
Then 


p(an) + p(all) ^ ^ 

er df 


I (l ■ [l-P^Wp^u,]" 2 ), 


Note that, for any combination of symbol errors introduced by the 
C^-decoder in the n^-l subsections of the u-th section of a data-block v, there 
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Is exactly one symbol error pattern introduced by the (^-decoder in the last 
subsection such that the decoded u-th section is a codeword in C^. The proba- 
bility that a given nonzero error pattern is introduced by the C^-decoder in a 
specific subsection of the u-th section for 1 5 u ^ m^ is upper bounded by the 


following expression: 


d -1 d -1-ij 

2 t. i 2 If 


ep 


0 7 5 


1 i"° 


i 2 -° 


n -d„ 
r 2 


KT 


i,+i. 


min{ Ld 2 -i 1 -i 2 -l)/2j ,n r -i 2 -d 2 ) 

£ 


h - 0 


n r' 1 2* d 2 


[p e (u >] [p 6 < u ’ 0) ] r 

\3 


n " d 2' i 2‘ h 


d 2'h 


max 

a*0 


L 

j-d 2 -i 1+ h-L(d 2 -i 1 -i 2 -l)/2j 


V 1 1 

. j . 


(22) 


It follows from (22) that 


m. 


,(all) 


er 


m 2 -l 


* 5 kr } f ten 1 

u-l ' i-1 ^ 1 J 


[l-pW<u) - P<f(u)] 


m 2 -i-l 


(23) 


Let p^ a11 ^ denote the right-hand side of (23). Then P^ 11 ^ serves as an upper 
er er 

bound on the probability of an incorrect decoding for a data block v. 


Overall Error Performance and Throughput Efficiency 

Again we assume that the code C r is very powerful so that at most one 

( 2 ) ( 2 ) 

retransmission is needed to recover a data-block. Let P &r and P^ be the 

probabilities of a decoding error and a decoding failure for the system respec 

( 2 ) ( 2 ) 

tively. Then we have the following bounds on P gr and P^ : 
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( 24 ) 


P (2) < P +P P +P (all) < P + (P +P )P +P (all) A P (2) 
er er es er er er v es er' er er - er 


P< 2) S P<‘ U > S P<* ll > + P<J U > A P< 2 > . 
df df er df - df 


(25) 


The system throughput is measured by 


k k 

(2) A K 1 K 2 


n l n 2 (1+P es ) 


(n 2 *d 2 +l)k 1 

n.n 9 (l+P +P ) 
1 Z er es 


(26) 


6 . Specific Schemes for NASA Near Earth Satellite Communications 

For NASA near earth satellite communications, we propose two specific 
schemes. For the first scheme, we choose n^-k^-£-8. No inner code is used. 

g 

The outer code C 2 is the extended (256,224) Reed-Solomon (RS) code over GF(2 ) 

(or a shortened version of this code). The (256,224) RS code is actually the 
NASA standard code for TDRS Systems with an additional information symbol. This 
code has 32 parity- check symbols and is capable of correcting any combination of 
t or fewer symbol errors and e or fewer symbol erasures provided that 2t+e < 32. 
Note that the length of this code, 256, is a multiple of 32. The code C^. is the 
shortened (64,32) RS code obtained from shortening C 2> is capable of cor- 

recting 16 symbol errors and is extremely powerful. Therefore, even in a very 
noisy situation, a transmitted data block should be recovered at most with one 
retransmission. 

The error performance of (^-decoder, C^-decoder and the overall system is 

given in Table 1 and Figures 4-8 for various values of t 2 (the designed error 

.2 

correcting capability of C 2 ) . For channel bit-error rate €-10 , we see from 

Table 1 and Figure 8 that the probability of a decoding failure of the system 

with at most one retransmission is upper bounded by 4x10*^. The probability of 

a decoding error (from Table 1 and Figure 7) of the system with t^-1 is upper 

bounded by 2x10 The system throughput is shown in Figure 9. For bit-error 

-2 

rate c-10 , the system throughput efficiency is about 50%. For channel bit- 

-3 

error rate €-10 , the probability of a system decoding failure is upper bounded 
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-20 

1.64 xlO ! The system throughput efficiency is simply the system code rate, 

-30 

95%, and the probability of a decoding error is upper bounded by 10 ! 

-3 

Therefore, for e < 10 , the system is practically error-free and the system 

throughput is nearly 100% of the overall code rate 

For the second specific scheme, we choose i-8 and m.j-6. The inner code 
is a distance-4 shortened (55,48) Hamming code which is used for correcting 
single error and detecting double errors in a frame. The base outer code C ^ is 

O 

again the (256,224) extended RS code over the Galois field GF(2 ). is 
interleaved by a depth of 6. The code is again the (64,32) shortened RS code 
over Gf (2®) . 

The error performance of the combined and C ^ decoders is given in 

Figures 10 and 11. The error performance of the decoder is shown in 

Figure 12. The overall error performance of the system is given in Figures 13 

and 14. Throughput efficiency of the system is shown in Figure 15. From these 

figures, we see that the second specific scheme also performs extremely well for 

-2 

channel bit error rate e < 10 

Since the second scheme uses an inner code and the base outer code is 

interleaved, it provides better performance than the first scheme. From Figures 

7 and 13, we see that the second scheme gives smaller error probability for the 

same channel bit-error rate t . Figures 16 and 17 gives a comparison of two 

schemes in probability of decoding failure and throughput efficiency. For 

. 2 

bit-error rate e less than 2 x 10 , the probability of a decoding failure of 

the second scheme is much smaller than that of the first scheme. As a result, 

the second scheme provides higher throughput efficiency than the first scheme 

-3 

(see Figure 17) for a certain range of bit -error rates. However, for £ < 10 , 

the first scheme gives higher throughput efficiency because, for this error 
rate, the probability of decoding failure of the first scheme is also extremely 
small and it uses less parity-check bits. The second scheme is more complicated 
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to implement and requires more buffer storage at both transmitter and receiver 
due to interleaving. 

7. Conclusion 

In this report, a robust error control coding scheme has been presented and 
analyzed. The scheme is a cascaded FEC scheme supported by parity retransmis- 
sions for further error corrections in parts (or subsections) . The extra 
parity-check symbols for further error correction are transmitted only when they 
are needed. When the channel is quiet or not so noisy, the scheme behaves like 
a conventional FEC scheme with throughput efficiency equal to the overall code 
rate. When the channel is noisy, parity retransmission is requested and extra 
parity symbols for correcting errors in subsections of a codeword are trans- 
mitted. The parity symbols sent in a retransmission are formed based on the 
data in the original transmission and a half- rate invertible code C r . These 
parity symbols contain the same amount of information as the original data. As 
a result, they can be used to recover the original data either by inversion or 
by decoding based on C^. If C f is powerful enough, at most one retransmission 
is needed to recover the erroneous data. If selective-repeat ARQ is used for 
retransmissions, the system throughput efficiency would remain high, say 50% of 

the system code rate, even for very high bit-error rate, say in the range of 
• 2 

10 . The proposed error control scheme is particularly suitable for satellite 

links with long propagation delay and nonstationary bit-error rate. The scheme 
uses the same amount of buffer store at the transmitter and receiver as a 
corresponding conventional hybrid ARQ scheme [2,3]. 

Two specific schemes are proposed for NASA's near earth satellite communi- 
cations. The first scheme does not use an inner code and the outer code is not 

g 

interleaved. The outer code is the NASA standard (256,224) RS code over GF(2 ). 
The half- rate code used for correcting errors in subsections during parity 
retransmission is the shortened (64,32) RS code obtained from the (256,224) RS 
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code. The (64,32) code is used for correcting 16 symbol errors over a subsec- 
tion pair of 64 symbols. Therefore, it is a very powerful code. As a result, 

an erroneous data would be recovered with at most one retransmission even for a 

-2 

bit-error rate as high as 10 . The scheme provides very high reliability (i.e. 

very low error probability) for channel bit-error rate up to 10 . If this 

scheme is used in cooperation with Yu-Lin's selective- repeat ARQ retransmission 

protocol, a buffer with size equal to the number of codewords transmitted in one 

round-trip delay is needed at the receiver. This buffer size is the same as the 

buffer used in a corresponding conventional selective -repeat ARQ scheme. 

The second specific scheme proposed for NASA uses the same (256,224) RS 

code as the base outer code C ^ except that it is interleaved to a depth of 6. 

The half -rate code C for parity retransmission is again the (64,32) shortened 

RS code. The scheme uses an inner code which is the distance-4 (55,48) 

shortened cyclic Hamming code. The inner code is used for correcting single 

error and detecting double errors in a frame. The decoding of this code is very 

simple, it can be implemented with a ROM using table-look-up. The second scheme 

provides better performance than the first scheme, however it requires a buffer 

with a size which is six times larger than that of the first scheme. For 

.3 

channel bit -error rate less than 5x10 , both schemes practically provide 

error- free communication and have the same throughput performance. However for 

-3 -2 

bit- error rate in the range 5x10 to 10 or higher, the second scheme gives 

much higher throughput efficiency. If NASA's satellite links for high rate file 

-3 

transfer operate with error rates less than 5x10 , we recommend that the first 

-2 

scheme be used. If 50% throughput efficiency at bit-error rate e-10 is 

acceptable, the first scheme is still a better choice. The second scheme is 

recommended when extremely high reliability and high throughput are needed for 

-3 -2 

bit-error rate in the range from 5x10 to 10 
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APPENDIX A 


Invertible Property of the Code 

Since C 2 is a maximum- distance -separable code, is also a maximum- 
distance -separable code [11]. The generator matrix in systematic form of C r is 
a (d 2 -l)x(2d 2 -2) matrix of the following form: 

G r - [I P] 

where I is a (d 2 *l) -dimensional identity matrix and P is a (d 2 -l)x(d 2 -l) square 

2 

matrix over GF(2 ). The parity-check matrix of C in systematic form is then 
T 

where P is the transpose of P. Since the minimum distance of is d 2 , the 
T 

d 2 *l columns of P must be linearly independent. This implies that the d 2 *l 
rows of P are linearly independent. As a result, P is a nonsingular square 

9 . l -1 

matrix over GF(2 ). Let P be the inverse of P. Then PP - I. 

_ 9 

Let a be a (d 2 *l) -tuple over GF(2 ). Then the (2d 2 *2) -tuple , 

v - (a,b) - a-G r , (A-l) 

is a codeword in C where 

r 

b - a-P (A-2) 


is the parity part of v. 

- & - - 
Let a' be another (d 2 *l) -tuple over GF(2 ) such that aVa. The codeword in 

C^_ corresponding to a' is 

v' - (i'.b') - I' -G r 

where 


b' - a' -P . (A-3) 

Now we want to show that b ** b'. Suppose that b-b'. Then, from (A-2) and 
(A-3), we have 

0 - b-b' - (a-a')P . (A-4) 


- 23 - 



Since a * a' , a-aV 0. Equation (A-4) implies that the rows of P are not 
linearly independent. This is a contradiction. Hence b ** b' . This says that 
there is a one-to-one correspondence between a and b. Multiplying both sides of 
(A- 2) by P \ we have 

b-P* 1 - a-P-P' 1 - 7 . 

Hence the data a can be recovered uniquely from its corresponding parity part b 


by taking an inversion. 



APPENDIX B 


Most commonly used maximum- distance -separable codes are RS codes. Let C, 


be a RS (or shortened or extended) code over GF(2 ). Let g 2 (X) be the generator 

polynomial of C 2 . The degree of is n 2~^2 “ d 2'^‘ Ijet *> e a code 

polynomial in C„. Let v 1 (X),v^(X) v (X) be the m„ subsections of v(X) . 

Z 1 Z m 2 l 


Then 


where 


_ r V k 9 _ (nu-l)(n 9 -k„) 

v(X) - v. (X) + v 0 (X)X 4 + ... + v (X)X * 

12 m 2 


V i (X) ’ V i0 + V il X + 


+ v. . .X 

i,n 2 -k 2 -l 


t»2 “^2 — T 


(B-l) 


(B-2) 


Since C r is obtained by shortening C^, the generator polynomial of C r is also 
g 2 <X). If we encode the i-th section v^(X) of v(X) based on C^, the parity 
check part of the code polynomial is the remainder R r (v^(X)] obtained by 

V k 2 - 

dividing X v^X) by g^X), i.e., 

V k 2 - - 

X Z v^X) - « 1 (X)g 2 (X) + R r [v i (X)] . (B-3) 

where the degree of R r (v^(X)] is n 2' k 2"^ or less. Consider the polynomial, 

_ _ _ r 

R[v(X)] - R r [v 1 (X)] + R r [v 2 (X)] X 


n 2' k 2 


+ R [v (X)]X 
r l m 2 


(® 2 *1) ( 1 * 2 ^ 2 ) 


(B-4) 


From (B-l), (B-3) and (B-4), we have 

- I n 2* k 2 

R[v(X)l - a^X) + a 2 (X)X z + ... 

(m 9 -l)(n 9 -k 9 )] _ _ n 9 -k 9 

+ a (X)X Z ) g 2 (X) + v(X)X Z Z . (B-5) 

Since v(X) is code polynomial in C^, v(X) is divisible by g^(.X) . From (B-5), we 
see that R[v(X)] is divisible by g 2 (X). Since R[v(X)] is a polynomial of degree 
m 2 (n 2 -k 2 )-l - ^-1 or less, R[v(X)] is a code polynomial in C 2 - 
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APPENDIX C 


Derivation of p^ max ^ 


-52- 


Let P^(w) be the probability of occurrence of an error pattern of weight w 


in a subsection after C^-decoding. It is shown in [5,6] that this probability 


is upper bounded by 


min(t r ,n r -w) 


P (w) < P (w) A 
er er - 


n -w 
r 


h-0 


w 


5 LI 

J j-w+h-t 


. w+h-d, 


r T r in 

2 . ~ r 

p 1-p 

r e *e 


n -w-h r 


max{p (a)} 
a*0 ® 


where n r - 2(d 2 -l), t f - L(d 2 *l)/2j and 


max(p (a) } - « (1- e) 


i-1 


a»*0 


Next we want to show that, if 


Hi OilOl 

1-e 


d 2 +3 


4(d 2 -l) ’ 


then 


P (d.) > P (w) 
er 2 er 


for d 2 < w < n r> From (C-l), we have 


min{t ,n -w-1} 
r r 


P (w+1) - 
er 


h-0 

I w+l+h-d 


n -w-1 

w 

r 

w+r 

r 

k * . 

5 

i-w+h-t 
j r 

j+i 


-j wrirn-Urt r 

P 1-P 

L e J e 


21, - 


max{p (a) } 


ai'O 


j+1 


(C-l) 


(C-2) 


(C-3) 


(C-4) 


By comparing each term in the right-hand sides of (C-l) and (C-4), we see that, 



(C-5) 



obtain the upper bound given by (11). 
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Table 1 Error performance and throughput efficiency 
of the first specific scheme 
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Byte Demultiplexer 























Original transmission 


1st retransmission 


2nd retransmission 


3rd retransmission 




Figure 3. Alternate data-parity retransmission strategy 
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Figure 5 The sum of probabilities of a decoding failure and a decoding 

error of the C2~decoder for a received data-block or parity-block. 
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Figure 6 Upper bounds on the probability of incorrect decoding for a 

block when each subsection pair is decoded by the C r -decoder. 
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, Figure 8 Upper bounds on the probability of a decoding failure of 
the first specific scheme. 
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Figure 9 Lower bounds on a measure of the throughput efficiency of 
the first specific scheme. 
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BIT-ERROR- RATE 

Figure 10 Upper bounds on the probability of a decoding error for a 
received data-block or parity-block by and C 2 _ decoders. 
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BIT-ERROR-RATE 

Figure 11 The sum of probabilities of a decoding failure and a decoding 
error for a received data-block or parity-block by the C, 
and C2~decoder . * 
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1 % 2 % 3 % 

BIT-ERROR-RATE 

Figure 12 Upper bounds on the probability of a decoding error for a block 


when each frame is decoded by C^-decoder at the first stage and 
then each received subsection pair is decoded by C r -decoder at 
the second stage. 
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Figure 14 Upper bound on the probability of a decoding failure for the 
second specific scheme. 
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Figure 17 Lower bounds on a measure of the throughput efficiency for 
the two schemes. 


